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Abstract 

Accurate calculations of the binding energies, the hyperfine split- 
ting, the bound-electron g-factor, and the parity nonconservation ef- 
fects in heavy few-electron ions are considered. The calculations in- 
clude the relativistic, quantum electro dynamic (QED), electron-correlation, 
and nuclear effects. The theoretical results are compared with avail- 
able experimental data. A special attention is focused on tests of QED 
in a strong Coulomb field. 

1 Introduction 

Accurate calculations of heavy few-electron ions must be performed in the 
framework of the rigorous quantum electrodynamic (QED) formalism. The 
basic methods of quantum electrodynamics were formulated to the beginning 
of 1930's, almost immediately after creation of quantum mechanics. This 
theory provided description of such low-order processes as emission and ab- 
sorption of photons and creation and annihilation of electron-positron pairs. 
However, application of this theory to some higher-order effects gave infinite 
results. This problem remained unsolved until the late 1940's when Lamb 
and Retherford discovered the 2s — 2pi/2 splitting (Lamb shift) in hydrogen. 
This discovery stimulated theorists to complete the creation of QED since it 
was believed that this splitting is of quantum electrodynamic origin. First 
evaluation of the Lamb shift was performed by Bethe who used Kramer's 



idea of the mass renormalization. A consequent QED formalism was devel- 
oped by Dyson, Feynman, Schwinger, and Tomonaga. They showed that all 
infinities can be removed from the theory by so-called renormalization pro- 
cedure. The basic idea of this procedure is the following. The electron mass 
and the electron charge, which originally occur in the theory, are not directly 
measurable quantities. All physical quantities calculated within QED be- 
come finite if they are expressed in terms of the physical electron mass and 
charge, parameters which can directly be measured in experiment. All cal- 
culations in QED are based on the perturbation theory in the fine structure 
constant a ~ 1/137.036. The individual terms of the perturbation series are 
conveniently represented by so-called Feynman diagrams. 

Before the beginning of 1970's investigations of QED effects in atomic 
systems were mainly restricted to low-Z atoms such as hydrogen or helium 
(here and below, Z is the nuclear charge number). In these systems, in 
addition to the small parameter «, there is another small parameter, which 
is aZ. For this reason, all calculations of low-Z atoms were based on the 
expansion in a and aZ. 

A great progress in experimental investigations of heavy few-electron ions, 
which was made for the last decades (see [1, 2] and references therein), has re- 
quired accurate QED calculations for these systems. Investigations of heavy 
few-electron ions play a special role in tests of quantum electrodynamics. 
This is due to two reasons. First, in contrast to low-Z atoms, the parameter 
aZ is not small and, therefore, the calculations must be performed without 
any expansion in aZ [3]. Second, in contrast to heavy neutral atoms, the 
electron-correlation effects can be calculated to a high accuracy using the 
pertubation theory in the parameter \jZ. For this reason, the QED effects 
are not masked by large electron-correlation effects, as it takes place in neu- 
tral atoms. This provides an excellent opportunity to test QED at strong 
electric fields. 

The calculations of high-Z few-electron ions are generally based on per- 
turbation theory. To zeroth order, one can consider that electrons interact 
only with the Coulomb field of the nucleus. The interelectronic-interaction 
and QED effects are accounted for by perturbation theory in the parame- 
ters XjZ and a, respectively. This leads to quantum electrodynamics in the 
Furry picture. To formulate the perturbation theory for calculations of the 
energy levels, transition and scattering amplitudes, it is convenient to use 
the two-time Green function method [4]. For very heavy ions the parameter 
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\jZ becomes comparable with a and, therefore, all the corrections may be 
classified by the parameter a only. In the present paper, we consider the 
current status of these calculations. 

Relativistic units {^ — c— \) are used in the paper. 



2 Binding energies in heavy few-electron ions 
2.1 H-like ions 

To calculate the binding energy in a hydrogenlike ion, we may start with the 
Dirac equation, 

(a • p + /?m + 1/c(r))^(r) = E^/'(r) , (1) 

where Vc(r) is the Coulomb potential induced by the nucleus. For the point- 
nucleus case, this equation leads to the binding energy 

E^^ -mc^^- ^ , mc\ 2 



where i' = n + \J{j + 1/2)^ — {ctZy — {j + 1/2), n is the principal quantum 
number, and j is the total angular momentum. To get the binding energy 
to a higher accuracy, we have to evaluate the quantum electrodynamic and 
nuclear corrections. 

The finite-nuclear-size correction can be obtained by solving the Dirac 
equation with the potential induced by an extended-charge nucleus and tak- 
ing the difference between the energies for the extended and the point nu- 
cleus model. This can be done either numerically (see, e.g., Ref. [5]) or 
analytically [6]. To a relative accuracy of ~ 0.2%, in the range Z=l-100 the 
finite-nuclear-size correction is given by the following approximate formulas 
[6] 



AE^, = ^-^[l + iaZrfUaZ)]{2^-^y'rnc\ (3) 

{OiZYv? — !,^ , ^,9^ , ^xn/„aZ R \27 2 
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where 7 = -i/l — (aZ)'^, 



fis{aZ) = 1.380 - 0.162Q;Z+1.612(aZ)^ 
f2s{aZ) = 1.508 + 0.215aZ + 1.332{aZf , 
f2p,/^{aZ) = 1.615 + 4:.319aZ -9.152{aZf + n.87{aZf , 

and i? is an effective radius of tlie nuclear cliarge distribution defined by 



i?={3(0[l-4(«^r(^^-^)]} • (5) 



Tfie corresponding formulas for other states can be found in Ref. [6]. We 
note that, in contrast to the nonrelativistic case where the corresponding 
correction is completely defined by the root-mean- square radius, in heavy ions 
the higher-order moments of the nuclear charge distribution are required to 
determine the nuclear-size correction on a 1% accuracy level. For instance, 
in case of Z—92 to gain the precision ~ 0.2% it is necessary to know the 
moment (r^). 

The next corrections one should take into account are the QED correc- 
tions of first order in a. They are determined by the self-energy (SE) and 
vacuum-polarization (VP) diagrams (Fig. la,b). The contribution of the 
self-energy diagram (Fig. la) combined with the corresponding mass coun- 
terterm is given by 

/oo r r 

xDf,^{uj, xi - X2)a''^a(x2) -5m j dx ^/'„(x)t/'„(x) . (6) 

Here ipai'^) is the Dirac-Coulomb wave function of the state under considera- 
tion, ip = ip'^^o, G{uj, xi, X2) is the Coulomb Green function, D^,j{uj, xi— X2) is 
the photon propagator, = (1, a), and a is a vector incorporating the Dirac 
matrices. The first evaluation of the SE correction for heavy ions was per- 
formed by Desiderio and Johnson [7] who employed the method suggested by 
Brown, Langer, and Schacfer [8] . Later, Mohr [9] developed another method 
which allowed him to perform a high precision evaluation of this correction 
in the range Z = 10—110. An alternative approach to this problem was 
worked out in Refs. [10, 11, 12]. The most accurate calculations of the SE 
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correction to all orders in aZ were performed by Mohr [13] and by Indelicato 
and Mohr [14] for the point nucleus case and by Mohr and Soff [15] for the 
case of extended nuclei. The highest accuracy for low-Z ions was achieved 
by Jentschura et al. [16]. 

The formal expression for the vacuum-polarization correction (Fig. lb) 
is given by 

A£; = ^ c^u; y c^xi 1 dx2 Vl(xi)^^-^^[TrG'(a;,X2,X2)]V'a(xi) . (7) 

This expression is ultraviolet divergent. It can be renormalized by dividing 
into two parts. The first part corresponds to the first nonzero term in the 
potential expansion of the Coulomb Green function in powers of aZ. This 
part, so-called Uehling part, becomes finite due to the charge renormalization 
and its evaluation causes no problem. The second part, so-called Wichmann- 
Kroll (WK) part, accounts for all higher-order terms of the «Z-expansion. 
Despite this part is finite, the regularization is still needed due to a spurious 
gauge-dependent piece of the light-by-light scattering contribution. Calcula- 
tions of the WK contribution in a wide range of Z were performed first by 
Soff and Mohr [17] for the extended nucleus case and by Manakov et al. [18] 
for the point nucleus case. The most accurate results for some specific ions 
were obtained by Persson et al. [19]. 

The QED corrections of second order in a are determined by diagrams 
presented in Fig. 2. Most of these diagrams can be calculated using the 
methods developed for the first-order SE and VP corrections [20, 21, 22, 23, 
24, 25, 26]. The most demanding problem is to evaluate the SE-SE diagrams 
[27] and the V(SE)P diagram (the diagram with a photon line inside the 
VP loop) [28]. The loop-after- loop SE-SE diagram was first evaluated by 
Mitrushenkov et al. [29] . A partial evaluation of the other SE-SE diagrams 
was performed by Mallampalli and Sapirstein [30] . The residual SE-SE terms 
were first calculated by Yerokhin and Shabaev [31]. Finally, in Refs. [32] the 
whole gauge invariant set of the SE-SE diagrams was evaluated in the range 
Z = 10 — 100. As to the S(VP)E diagram, to date it was evaluated only in 
the Uehling approximation [20, 21]. 

The calculations discussed above are based on the approximation in which 
the nucleus is considered as a source of the external Coulomb field. It defines 
quantum electrodynamics within the external field approximation. The first 
step beyond this approximation would consist in evaluation of the nuclear 
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recoil correction. As is known, in the nonrelativisic theory of a hydrogenhke 
atom the recoil effect can easily be taken into account by using the electron 
reduced mass — mM/{m + M). A rigorous relativistic theory of the re- 
coil effect can be formulated only in the framework of QED. In case of a 
hydrogenlike atom, a close formula for the recoil effect to first order inm/M 
and to all orders in aZ was derived in Rcf. [33] (sec also [34] and references 
therein). According to this formula, the recoil correction is given by the sum 
of a low-order term AE'l and a higher order term A£^h, where 

A^^L = ^(a[[p2-(D(0)-p + p-D(0))][a), (8) 

xG(u; + £;„)(DH + |^)|a). (9) 

Here p is the momentum operator, G{uj) is the Coulomb Green function, 
Dmi!^) — ~4='n'aZaiDim{uj) , and Dik{u!,r) is the transverse part of the pho- 
ton propagator in the Coulomb gauge. In equation (9), the scalar product 
is implicit. The term AE'l contains all the recoil corrections within the 
(aZ)^m^/M approximation. For the point nucleus case, it can easily be 
calculated analytically [33] 

A^, = !^4^. (10) 

The term AEn contains the contribution of order {aZ)^m'^ /M and all con- 
tributions of higher order in aZ which are not accounted for by the AEl 
term. Numerical evaluations of the recoil correction to all orders in aZ were 
performed in Refs. [35, 36] for point and extended nuclei, respectively. 

Finally, one should take into account the nuclear polarization correction, 
which sets the ultimate limit of accuracy up to which QED can be tested in 
atomic systems. This correction is determined by the electron-nucleus inter- 
action diagrams in which the intermediate states of the nucleus are excited. 
It was evaluated by Plunien and Soff [37] and by Nefiodov et al. [38]. 

The individual contributions to the ground-state binding energy in 

238U91+ 

are given in Table 1. The uncertainty of the Dirac binding energy comes 
from the uncertainty of the Roohc constant [39]. The finite-nuclear-size cor- 
rection is evaluated for the Fermi model of the nuclear charge distribution 
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with (r^)^/^ = 5.8507(72) fm [40]. The uncertainty of this correction is esti- 
mated by adding quadratically two errors, one obtained by varying the root- 
mean-square radius and the other obtained by changing the model of the 
nuclear-charge distribution from the Fermi to the homogeneously- charged- 
sphere model. As one can see from the table, the present status of the theory 
and experiment on heavy H-like ions provides a test of QED on the level of 
about 2%. 

2.2 Li-like ions 

To date, the highest accuracy was achieved in experiments on the 2pi/2,3/2— 2s 
transitions in heavy Li-hke ions [42, 43, 44, 45]. In these systems, in addi- 
tion to the one-electron contributions discussed above, one has to evaluate 
two- and thrcc-clcctron contributions. To first order in a, the two-clcctron 
contribution is determined by the one-photon exchange diagram (Fig. 3) 
whose calculation causes no problem. To second order in a, one should ac- 
count for the two-photon exchange diagrams (Fig. 4) and the self-energy 
and vacuum-polarization screening diagrams (Fig. 5). The two-photon ex- 
change diagrams were evaluated by different authors [46, 47, 48, 49, 50]. 
The results of these calculations are in a good agreement with each other. 
The complete calculation of the vacuum-polarization screening diagrams was 
performed in Ref. [51]. The self-energy screening diagrams were evaluated 
in Refs. [52, 47]. To gain an accuracy required by the experiments, in 
addition to higher-order one-electron QED corrections, one should evaluate 
the interelectronic-interaction corrections of third and higher orders in the 
parameter \ jZ. Such cvahiations within the framework of the Breit approx- 
imation were accomplished in Refs. [53, 54, 48]. 

The individual contributions to the 2pi/2 — 2s transition energy in Li- 
like uranium are presented in Table 2. The total theoretical value of the 
transition energy, 280.66(22) eV, agrees well with the related experimental 
values, 280.59(10) eV [42] and 280.52(10) cV [45]. Comparing the first- and 
second-order QED contributions, —42.93 eV and 1.45(18) eV, respectively, 
with the total theoretical and experimental uncertainties, we conclude that 
the present status of the theory for Li-like uranium provides a test of QED 
on a 0.5% level to first order in a and on a 15% level to second order in a. 
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2.3 He-like ions 



In Refs. [55, 56] the two-electron contribution to the ground state energy 
of hehumUke ions was measured. It was done by comparing the ionization 
energies of heliumlike and hydrogenlike ions. This experiment is of special 
importance since to date the two-electron contribution is the only measured 
value which has been calculated to the second order in a. As in case of 
Li-like ions, the lowest-order two-electron contribution is determined by the 
one-photon exchange diagram. The contributions of second order in a are 
given by the two-photon exchange diagrams and by the SE and VP screen- 
ing diagrams. The two-photon exchange contribution was first evaluated by 
Blundell et al. [57] and by Lindgren et al. [58]. This contribution can be 
conventionally divided into two parts: one which corresponds to the Breit 
approximation and the other which is beyond the Breit approximation. The 
VP and SE screening diagrams for the ground state were evaluated in Refs. 
[59, 60, 61, 62]. The corresponding calculations for excited states of He-hke 
ions were performed in Refs. [63, 64, 65, 48, 49, 66, 67]. 

In Table 3 we present the individual contributions to the two-electron 
binding energy of the ground state in He-like uranium. The two-photon 
exchange contribution is divided into two parts as described above. The 
Breit contribution and the higher-order QED corrections are evaluated as in 
Ref. [67]. As one can see from the table, to test the screened QED effects, 
the experimental precision should be improved by an order of magnitude. 

3 Hyperfine splitting in heavy ions 

High-precision measurements of the hyperfine splitting (HFS) in heavy hy- 
drogenlike ions [68, 69, 70, 71, 72] have triggered a great interest to theo- 
retical calculations of this effect. The ground-state hyperfine sphtting of a 
hydrogenlike ion is conveniently written as [73] 



Here is the proton mass, n is the nuclear magnetic moment, /x^v is the 
nuclear magneton, and / is the nuclear spin. A{aZ) denotes the relativistic 



AE^ = -a{aZ) 




(11) 
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factor 



A{aZ) 



1 



l + ^-(aZr+^-l(aZr + ..., 



(12) 



7(27 - 1) 



S is the nuclear charge distribution correction, e is the nuclear magnetization 
distribution correction (so-called Bohr-Weisskopf correction), and Xj-ad is the 
QED correction. The most accurate calculations of the QED corrections were 
performed in Refs. [74, 75, 76]. The uncertainty of the theoretical predictions 
is mainly determined by the uncertainty of the Bohr-Weisskopf (BW) effect. 
In calculations, based on the single-particle nuclear model [73, 77, 74, 78], 
which provide a reasonable agreement with the experiments (see, e.g., Ref. 
[4]), this uncertainty may amount up to about 100% of the BW effect and 
is generally larger than the total QED contribution. More elaborated cal- 
culations, based on many-particle nuclear models [79, 80], do not provide a 
desirable agreement with the experiments. 

A new method to determine the BW effect was recently developed in Ref. 
[81]. In this method the BW correction to the hyperfine splitting in hydro- 
genhke ^°^Bi, ^osrp^^ 205'p| determined using experimental data on 
the hyperfine splitting in the corresponding muonic atoms [82, 83]. The pa- 
rameters of the nuclear magnetization distribution were chosen to reproduce 
the experimental values of the nuclear magnetic moment as well as the BW 
effect in muonic atoms extracted from the corresponding experiments. The 
single-particle and configuration-mixing nuclear models were considered. To 
increase the precision of determining the BW contribution, the QED correc- 
tions for muonic atoms have been evaluated. In Table 4, the BW correc- 
tion to the HPS in electronic H-like ions, derived from the experiments on 
muonic atoms, is compared with the BW correction obtained by direct cal- 
culations within the single- and many-particle nuclear models. Taking into 
account that the uncertainty of the single-particle results may amount up to 
30-100%, we conclude that the e values, based on experiments with muonic 
atoms, have a better accuracy. 

In Table 5, we compare the total theoretical results for the HFS in H-like 
ions with experiment. As one can see from the table, the results based on ex- 
periments with muonic atoms are closer to the experimental ones, compared 
to the results based on the direct calculations within the single-particle nu- 
clear model. However, due to a higher accuracy of the former results, a small 
discrepancy between the theory and experiment occurs for 209gj82+ 
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also be seen that in case of two isotopes of Tl the results of Ref. [81] are much 
closer to the related experimental values than in case of Bi. Since, in contrast 
to the Tl isotopes, the bismuth nucleus has a nonzero electric-quadrupole 
moment, one may expect that the discrepancy between the theory and ex- 
periment is caused by a large contribution of the electric-quadrupole splitting 
in muonic bismuth. It is known that the first-order electric-quadrupole split- 
ting vanishes for s states. However, the second-order effect is nonzero and, 
in principle, it may be significant for muonic atoms. This is due to a rela- 
tively large role of the electric-quadrupole HFS interaction in muonic atoms 
compared to electronic ones. Our numerical evaluation of this effect showed 
that, whereas it changes the HFS in muonic bismuth by about —3.4 eV, it 
almost does not affect the BW value presented in Table 4. 

One of the main goals of the HFS experiments with heavy H-hke ions was 
to probe the magnetic sector of QED in the presence of a strong Coulomb 
field. The analysis of the theoretical results showed, however, that the QED 
corrections to the HFS in heavy H-like ions are strongly masked by the un- 
certainty of the BW effect. This makes unfeasible tests of QED by a direct 
comparison of theory and experiment on the HFS in heavy H-like ions. An 
opportunity for QED tests has been found by considering a specific difference 
of the ground-state HFS values in H- and Li- like ions [84]. Namely, it was 
shown that the difference 

A'E = AE^^'^ - CAE^^'^ , (13) 

where AE^^'^ and AE^^'^ are the HFS in H- and Li-like ions of the same 
isotope, is very stable with respect to variations of the nuclear model, if 
the parameter ^ is chosen to cancel the BW corrections in the right-hand 
side of equation (13). The parameter ^ is almost independent of the nuclear 
structure and, therefore, can be calculated to a high accuracy. In case of 
Bi, the calculations yield ^ = 0.16885 and A'E= -61.27(4). The non-QED 
and QED contributions amount to -61.52(4) and 0.24(1), respectively. It 
means that the QED contribution is six times larger than the current total 
theoretical uncertainty. This provides good perspectives for tests of QED in 
the HFS experiments. In particular, it can be shown that this method allows 
one to test QED on level of a few percent, provided the HFS is measured to 
accuracy ~ 10~^. 
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4 Bound-electron g-factor 



The g-factor of a hydrogenlike ion with a spinless nucleus can be defined as 

_ iJM,mJM.) ^ (14) 

IJ,bMj 

where /x*^*^^ is the operator of the magnetic moment of electron and is the 
Bohr magneton. For the Is state, a simple relativistic calculation based on 
the Dirac equation yields 

^D = 2-(4/3)(l-^l-(aZ)2). (15) 
The total g-factor value can be written as 

g^'^ = + A^QED + A^rec + A^NS , (16) 

where Aqqed is the QED correction, Af^rec is the nuclear recoil correction, 
and Ag^NS is the finite-nuclear-size correction. High-precision measurements 
of the bound-electron g-factor in H-like carbon [85, 86] and oxygen [87] have 
stimulated successful theoretical calculations of this effect (see Refs. [88, 
89, 90, 91, 92, 93, 94] and references therein). In particular, these studies 
provided a new determination of the electron mass to an accuracy which 
is four times better than that of the 1998 CODATA value. As a result, 
the 2002 CODATA value for the electron mass [39] is derived from the g 
factor measurements. An extension of these experiments to higher- Z systems, 
which is anticipated in the near future, could lead also to an independent 
determination of the fine structure constant. 

For heavy H-like ions the theoretical uncertainty of the bound-electron 
g-factor is mainly determined by the nuclear effects [91, 93]. This uncertainty 
becomes comparable with the binding QED corrections of second order in a 
and, therefore, strongly restricts probing QED in these investigations. How- 
ever, in Ref. [95] it was shown that the uncertainty caused by the nuclear 
effects can be significantly reduced in a specific difference of the g-factors of 
H- and Li-hke ions of the same isotope. This gives a good opportunity for 
tests of the magnetic sector of QED in the presence of a strong Coulomb 
field. The most accurate calculations of the g-factor of Li-like ions in the 
range Z=6— 92 were performed in Ref. [96]. 
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The g-factor of a H-like ion with nonzero nuclear spin can be approxi- 
mated as 



^ ^ j,) F{F+l)+j{j + l)-I{I + l) 
^ ^ 2F(F + 1) 

m F(F+l) + /(/ + l)-j(j + l) 

2F(F + 1) ■ ^^^^ 

Here m and rUp are the electron and the proton mass, respectively, g'^'^^ is the 
bound-electron g-factor defined above, gi = n/Hnl is the nuclear g-factor, 
A*n = \e\/2mpC is the nuclear magneton, j and / are the electron and the 
nucleus angular momentum, respectively, and F is the total angular momen- 
tum of the ion. According to formula (17) the contribution of the nuclear 
g-factor is suppressed by factor m/mp compared to the electronic g-factor. It 
follows that measurements of the g-factor of ions with nonzero nuclear spin 
with a 10~^ accuracy would provide determinations of the nuclear magnetic 
moments on the 10~^ accuracy level. Calculations of various corrections to 
formula (17) were presented in Ref. [97]. 

Another possibility for investigations of the g-factor of ions with nonzero 
nuclear spin was discussed in Ref. [98]. In that work it was shown that the 
transition probability between the ground state hyperfine splitting compo- 
nents of a hydrogenlike ion, including the QED and nuclear corrections, is 
given by 



3 



w — — 



3 21 + 1 



m 



r ' - 9i — 

nip 



(18) 



where uj is the transition frequency. Formula (18) allows one to calculate 
the QED and nuclear corrections to the transition probability using the cor- 
responding corrections to the bound-electron g-factor. In particular, it was 
found [98] that in cases of Pb and Bi the QED and nuclear corrections in- 
crease the transition probability by about 0.3%. In Ref. [99] the life time 
of the upper hyperfine splitting component in 209gj82+ measured to be 
''■exp=397.5(1.5) fxs. This result is in a good agreement with the theoretical 
prediction [98] Ttheo = 399.01(19)//s. Using formula (18) and the experimen- 
tal values of the hyperfine splitting and the transition probability in 209gj82+ 
[68, 99], one finds the experimental value of the bound-electron g-factor in Bi 
to be 1.7343(33). The corresponding theoretical value is 1.7310. The values 
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of the individual contributions to the bound-electron g-factor in Bi are given 
in Table 6. From this table, it can be seen that including the QED correction 
is needed to obtain the agreement between the theory and experiment. 

5 Parity nonconservation effects with heavy 
ions 

Investigations of parity nonconservation (PNC) effects in atomic systems play 
a prominent role in tests of the Standard Model (SM) [100]. The well-known 
cesium experiment by Wieman's group [101], compared to the most elabo- 
rated theoretical result (see Ref. [102] and references therein), provided the 
most accurate test of electroweak theory at the low-energy regime. Further 
improvement of tests of the Standard Model with neutral atoms, from theo- 
retical side, is mainly limited by difficulties with accurate calculations of the 
electron-correlation effects. For this reason, the PNC experiments with heavy 
few-electron ions, where the electron-correlation effects can be evaluated to a 
high accuracy, seem highly desirable. The PNC effects in these systems were 
first discussed in Refs. [103, 104, 105]. For tests of the spin-independent part 
of the week interaction, a promising situation occurs in heavy He-like ions 
with the nuclear charge number near Z = 64 and Z = 92, where two levels of 
the opposite parity, 2^5*0 and 2^Pq, are very close to each other. The study 
of PNC effects with these ions requires precise knowledge of the 2^5*0 — 2^Po 
energy difference. The most accurate calculations of this difference, which 
include a complete set of two-electron QED corrections, were performed in 
Ref. [67] . In Table 7, the results of Ref. [67] are compared with the related 
data by other authors. 

A feasible PNC experiment with heavy ions was suggested by Labzowsky 
et al. [108]. Instead of the standard measurement of the circular dichroism 
which is rather difficult to perform with x-ray radiation, it was proposed 
to study a quenching-type experiment with interference of hyperfine- and 
weak-quenched transitions in polarized He-like europium. The idea of the 
experiment is the following. The basic one- photon decay 2^5*0 — 1^5*0 channel 
is the hyperfine-quenching (hfq) Ml transition which is due to the hyperfine 
mixing of the 2^5*0 and 2^5*1 levels. Another one-photon decay 2^5*0 — 1^5*0 
channel is due to the mixing of the 2^5*0 and 2^Pq levels caused by the weak 
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interaction of electrons with the nucleus and due to the hyperfine mixing of 
the 2^Po and 2^ Pi levels. As a result, the total amplitude of the one-photon 
2^5*0 — l^So transition is a mixture of the basic Ml and the additional El 
amplitude. The polarization of Eu^^"*" in the 2-'^5'o state can be described by 
the density matrix [108] 



1 



P 



3A 



1 + 7^(^-1) > (19) 



27 + 1 L 7 + 1 



where I is the operator of the total angular momentum, u is the unit vector 
directed along the ion polarization, and Aq is the degree of polarization (Aq < 
1). The probability for the emission of a photon in direction n is given by 
[108] 

dW{n)^^[l + e{u-n)]dQ, (20) 

where Wmi is the total probability of the hfq Ml 2^5*0 — 1^5*0 transition and 
e is the asymmetry coefficient caused by the PNC effects. Our evaluation of 
this coefficient yields 

e = QXo^W^i/Wui/il + 1) , (21) 

where Wei is the total probability of the weak-hyperfine quenching El 2^So — 
I^Sq transition. We note that expression (21) differs by about a factor of 2 
from the related expression derived in Ref. [108]. 

The experiment should consist in observing the difference of the transi- 
tion probability (20) due to a change of the ion polarization direction or, 
equivalently, due to rotating the detector around the beam direction by an 
angle tt. The value of this difference is proportional to the asymmetry pa- 
rameter e. Our calculation of this parameter for europium employing the 
transition energies from Ref. [67] (sec Table 7) yields e ~ 0.0004Ao. This 
value is almost four times larger compared to that obtained in Ref. [108]. 
This is due to a change of the 2^So — 2^Po energy difference, compared to that 
used in Ref. [108], and due to the change of the expression for e discussed 
above. 

The experiment under consideration requires preparing and storing a po- 
larized ion beam. A promising idea for preparation of polarized ion beams 
was suggested in Ref. [109]. 
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Conclusion 



In this paper we have considered the present status of calculations of the 
QED effects in heavy few-electron ions. Calculations of the PNC effects for 
heavy ions have also been discussed. 

To date, the most accurate tests of QED effects on binding energies in 
a strong Coulomb field have been accomplished in heavy Li-like ions: on 
a 0.5% accuracy level to first order in a and on a 15% accuracy level to 
second order in a. An improvement of the experimental accuracy by an 
order of magnitude accompanied by more accurate determinations of higher- 
order QED and nuclear-size effects would provide tests of QED beyond the 
external field approximation. 

High-precision measurements of the hyperfine splitting in heavy H-and 
Li-like ions of the same isotope are highly desirable. They would give a 
unique opportunity for tests of the magnetic sector of QED in the strongest 
electric field currently available for experimental study. 

The QED theory of the bound-electron g-factor has been probed by di- 
rect measurements on hydrogenlike carbon and oxygen. These measurements 
have also provided a new determination of the electron mass with an accuracy 
which is four times better than that of the previously accepted value. Exten- 
tions of these measurements to higher-Z systems and to ions with nonzero 
nuclear spin would provide independent determinations of the fine structure 
constant and the nuclear magnetic moments. 

Investigations of the PNC effects with heavy ions seem very promising 
for tests of the Standard Model. 
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Table 1: Individual contributions to the ground-state binding energy in 
238^91+^ in eV. The Lamb shift is defined as a part of the binding energy 
that is beyond its point-nucleus value given by equation (2). 



Point-nucleus binding energy 


-132279.93(1) 


Finite nuclear size 


198.33(52) 


First order SE 


355.05 


First order VP 


-88.60 


Second-order QED 


-1.26(33) 


Nuclear recoil 


0.46 


Nuclear polarization 


-0.20(10) 


Lamb shift theory 


463.78(62) 


Lamb shift experiment [41] 


460.2(4.6) 



Figure 1: First-order self-energy and vacuum-polarization diagrams. 

(a) (b) 
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Table 2: The 2pi/2 - 2s transition energy in 238u89+^ 



One-photon exchange 


368.83 


One-electron nuclear size 


-33.27(8) 


First-order QED 


-42.93 


Two-photon exchange within 




the Breit approximation 


-13.54 


Two-photon exchange beyond 




the Breit approximation 


0.17 


SE and VP screening 


1.16 


Three- and more photon exchange 


0.16(7) 


Nuclear recoil 


-0.07 


Nuclear polarization 


0.03(1) 


One-electron second-order QED 


0.12(18) 


Total theory 


280.66(21) 


Experiment [42] 


280.59(10) 


Experiment [45] 


280.52(10) 



Table 3: The two-electron binding energy of the ground state in 238^90+^ 
eV. 



One-photon exchange 


2265.90(1) 


Two- and more photon exchange 




within the Breit approximation 


-11.96 


Two-photon exchange 




beyond the Breit approximation 


-0.85 


SE screening 


-9.78 


VP screening 


2.63 


Higher-order QED 


-0.05(18) 


Total theory 


2245.89(18) 


Experiment [56] 


2248(9) 
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Table 4: The Bohr-Weisskopf correction £, derived from experiments on 
muonic atoms [81], is compared with previous direct evaluations of this ef- 
fect within the single-particle [77, 78] and the many-particle [79, 80] nuclear 
model. 



209gj82+ 205'jj80+ 203'jj80+ 



Elizarov et al. [81] 
Shabaev [78] 
Labzowsky et al. [77] 
Tomaselli et al. [79] 
Sen'kov and Dmitriev 



0.0123(15) 
0.0118 
0.0131 
0.0210 
[80] 0.0095(+7, -38) 



0.0193(27) 0.0155(40) 
0.0179 0.0179 



Table 5: Hyperfine sphtting in H-like ions, in eV. 



Theory [81] Theory [78] Experiment 

203Tpo+ 3.220(20) 3.229(17) 3.21351(25) [72] 

3.238(9) 3.261(18) 3.24410(29) [72] 

209^^82+ 5.098(7) 5.101(27) 5.0840(8) [68] 



Table 6: The bound-electron g-factor in ^°^Bi' 



Point-nucleus Dirac value 


1.7276 


QED 


0.0029 


Nuclear size correction 


0.0005 


Total theory 


1.7310 


Experiment [99] 


1.7343(33) 
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Table 7: The 2^Pq-2^Sq transition energy, in eV. 



z 


Artemyev et al. 


Andreev et al. 


Plant e et al. 


Drake 




[67] 


[49] 


[106] 


[107] 


63 


-0.224(69) 


-0.591 




-0.168 


64 


0.004(74) 


-0.389 


-0.170 


0.067 


65 


0.32(12) 


-0.153 




0.328 


66 


0.495(84) 


0.016 


0.341 


0.614 


89 


1.61(46) 






1.731 


90 


0.61(46) 




-0.095 


0.718 


91 


-0.31(55) 


-1.971 




-0.209 


92 


-2.64(28) 


-4.511 


-2.639 


-1.816 



Figure 2: Second-order one-electron Feynman diagrams. 
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Figure 3: One- photon exchange diagram. 




Figure 4: Two-photon exchange diagrams. 




ure 5: Self-energy and vacuum-polarization screening diagrams. 
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